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Most biomolecules occur in mirror, or chiral, images of each other. However, life is homochiral: 
proteins contain almost exclusively L-amino acids, while only D-sugars appear in RNA and DNA. 
The mechanism behind this fundamental asymmetry of life remains an open problem. Coupling 
the spatiotemporal evolution of a general autocatalytic polymerization reaction network to external 
environmental effects, we show through a detailed statistical analysis that high intensity and long 
duration events may drive achiral initial conditions towards chirality. We argue that life's homochi- 
rality resulted from sequential chiral symmetry breaking triggered by environmental events, thus 
extending the theory of punctuated equilibrium to the prebiotic realm. Applying our arguments 
to other potentially life-bearing planetary platforms, we predict that a statistically representative 
sampling will be racemic on average. 



I. INTRODUCTION 

Louis Pasteur was the first to realize, in the late 1840s, that many biomolecules display a mirror asymmetry known 
as chirality [24]. Although solutions synthesized in the laboratory are always racemic, with equal concentrations of 
each chiral isomer (enantiomer), the same is not true of living matter: proteins contain almost exclusively L-amino 
acids, while only D-sugars appear in RNA and DNA [8 . Several mechanisms have been proposed to explain life's 
homochirality, although none is conclusive [M] . Chiral biasing effects from parity violation (PV) in the weak interac- 
tions of particle physics [U [20] are not viable [16], while prebiotic biasing from circularly polarized UV light (CPL) 
originating in active star- forming regions [T1[T6] remains controversial [6]. Here, we propose that the homochirality of 
life may be explained by extending to prebiotic times the punctuated equilibrium hypothesis of Eldredge and Gould, 
whereby speciation occurred through alternating periods of stasis and intense activity prompted by external influ- 
ences [12] . We note that we are borrowing the concept of punctuated equilibrium with some freedom: the network of 
chemical reactions we describe is a non-equilibrium open system capable of exchanging energy with the environment. 
The periods of stasis that develop correspond to steady-states in that even though environmental influences may be 
negligible, chemical reactions are always occurring so as to keep the average concentrations at a constant value. 

A much-debated question is whether the observed homochirality of biomolecules is a prerequisite for life's emergence 
or if it developed as its consequence [4 , 9 . It has been argued that achiral peptide nucleic acids can form base-pairs and 
helical structures and thus are attractive candidates for the first carriers of genetic information ^3] • Considering that 
very specific conformations of structural entities such as a-helices and /3-sheets can only form from enantiomerically 
pure building blocks [8] [14] , we here take the bottom-up approach and assume that the conditions for the subsequent 
development of complex biomolecules had to be chiral. A related question is whether chiral compounds were fed into 
the prebiotic soup from outer space \J_^. The discovery of an excess of chiral organic compounds in the Murchison 
meteorite has supported this view [10l[TT]. Examination of isotopic distributions (in ^^N/^^N) and of a-branched 
amino acids of extraterrestrial origin has eliminated the possibility of contamination by Earth's biosphere, with 
measurements showing an excess of L-alanine of 50% and of 30% for L-glutamic acid [13 . A plausible explanation is 
that the cloud that originated the Solar System was subjected to CPL, such as synchrotron radiation from a neutron 
star. If this is indeed a viable mechanism, the same chiral bias should be prevalent throughout the Solar System 
but not necessarily throughout the galaxy. This should be contrasted with bias from PV, which should be the same 
throughout the universe. Finding an excess in D-chiral compounds in Mars [19j or elsewhere in the Solar System 
would contradict both scenarios. 

We propose an alternative point of view, that environmental effects might destroy any memory of a prior chiral 
bias, whatever its origin. Life's chirality is interwoven with early-Earths environmental history; specifically, with how 
the environment influenced the prebiotic soup that led to first life. 
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II. ENVIRONMENTAL EFFECTS ON SPATIOTEMPORAL POLYMERIZATION 



A. Polymerization Model 

We start with Sandars' polymerization model a generalization of Frank's pioneering approach [T5^, featuring 
autocatalysis with enantiomeric cross-inhibition. Consider a left-handed polymer L^, made of n left-handed monomers, 
Li. It may grow by adding another left-handed monomer with a rate kg, or be inhibited by adding a right-handed 
monomer Di with a rate kj. (Note that we denote D-compounds by the letter "D" as opposed to the notation set in 
Sandars' work where such molecules were denoted as "R".) The reaction network for n = 1, . . . , A/", where N is the 
maximum polymer length in the system, can be written as: 

r t 



Di + L„Di '^DiLnDu 



(1) 



supplemented by reactions for D-polymers by interchanging L ^ and by the production rate of monomers from 

the substrate: S Li; S Di. Cl(^d^ determine the enzymatic enhancement of L(I))-handed monomers, 

usually assumed to depend on the largest polymer in the reactor pool, Cl(^d^ = Ln{Dn) [25 , or on a sum of all 
polymers [28 . Soai's group obtained the best-known illustration of this autocatalytic mechanism with enantiomeric 
cross-inhibition [26 , with dimers {N = 2) as catalysts [3 . 

A set of coupled, nonlinear ordinary differential equations for the various concentrations, [Li], [l^i], [Ln], 
[Dn], describes the time evolution of the reaction network of eq. [l] supplemented by the equation for the substrate, 
d[S]/dt = Q{Ql -\- Qd)^ where Q is the substrate's production rate, and Ql — Qd = kcf[S]{CL — Cd) gives the net 
chiral excess in monomer production. / is the enzymatic fidelity, usually set to unity to maximize chiral separation. 
In Fig. [l] we show numerical solutions for polymerization reactions with A/" = 2, 5, and oo: starting as racemates, 
they all evolve toward homochirality. Also included is the solution for = 2 within the adiabatic approximation, 
where the rate of change for dimers and the substrate is assumed to be much slower than that of monomers, that is, 
when ksj « kc [5|. In Gleiser and Walker [18 , a detailed study of the polymerization reaction network for various 
values of N has shown that the trends displayed in Fig. [l] are true even when the effects of spatial dynamics are 
considered, as we will do next. Gleiser and Walker also concluded that although the adiabatic approximation predicts 
faster approach to steady-state conditions when compared with the full N = 2 model, it does produce the correct 
asymptotic values for the various concentrations. 



B. Spatiotemporal Polymerization 

We now extend this treatment to include spatial dependence, establishing a reaction-diffusion network by substitut- 
ing d/dt d/dt — kV'^^ where k is the diffusion constant [5]. In this coarse-grained approach, the number of molecules 
per unit volume is large enough so that the concentrations vary smoothly in space and time. The spatiotemporal evo- 
lution of the network is obtained by solving the coupled system of nonlinear PDEs for arbitrary values of n. Clearly, 
as n increases, solving and statistically analyzing the coupled system of equations in two and three spatial dimensions 
for various parameters (see below) becomes highly CPU intensive, even for the NCSA supercomputers we used. 

Given that Soai's reaction [26 - the only illustration so far of an autocatalytic network leading to chiral purity - 
featured dimers as catalysts, we focus here on the truncated system for A/" = 2 within the adiabatic approximation 
since, as with spatially-independent reaction networks (see Fig. [T]), it has similar qualitative behavior to networks 
with longer (A^ > 2) polymer chains [18]. The system then reduces to two coupled PDEs for the concentrations [Li] 
and [Di]^ being thus more amenable to a detailed statistical study while maintaining the key qualitative features of 
a larger reaction network. 

Introducing the variables S = X ^ Y and A = X — F, we can rewrite the reaction network in terms of the 
net dimensionless chiral asymmetry A{t^x^y^ z). For near-racemic initial conditions (|^(0, x, z)| < 10~^), the 
spatiotemporal evolution leads to the formation of left and right-handed percolating chiral domains separated by 
domain walls, as is well-known from systems in the Ising universality class (see Fig. [2|. Surface tension drives the 
walls until their average curvature matches approximately the linear dimension of their confining volume. At this 
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FIG. 1: The time evolution of the net chirahty 0{t) = ~ Dn)/ ^{Ln + Dn) from the reaction network of eq. [T]for several 

maximum polymer lengths [N). All models display similar qualitative behaviour, achieving pure chirality in short time-scales. 
The curves for N = 5 and N = oo overlap, lo = 



point, wall motion becomes quite slow, d{A{t))/dt 0, where {A{t)) is the spatially-averaged value of the net chiral 
asymmetry, and the domains coexist in near dynamical equilibrium in that the net stresses add to zero (see Fig. [2j 
top right). The time evolution of A{t) is shown in Figjsj For such model systems, it has been shown that the presence 
of a bias from PV or most GPL sources (even in the unlikely situation where they could be sustained unperturbed 
for hundreds of millions of years), would not lead to chirally-pure prebiotic conditions [16]. 



C. Coupling Environmental Effects 

Motivated by the experimental demonstration that stirring [2TJ [27j can effectively bias chirality, and a recent 
numerical study for a similar N = 2 system showing that the addition of fluid turbulence can speed up chiral evolution 
[5 , we proceed to examine the impact of random environmental effects in the reaction network. This is achieved in the 
simplest possible way via a generalized spatiotemporal Langevin equation [17]. The dynamical equations are written 
as: 



^o"' - kV^Sj = 1 - 52 + yj^t, x), (2) 

where Iq = {2ksQY^'^ ^ il'(x, t) is a dimensionless Gaussian white noise with two-point correlation function 
{w{-x' ,t')w{^,t)) = a^S{t' - t)^(x' - x), where is a measure of the environmental influence's strength. An Ising 
phase diagram can be constructed showing that (A) ^ for a > ac'. chiral symmetry is restored [17 . The value 
of ttc has been obtained numerically in two (a^ = 1.15(/c//o)cm^s) and three (a^ = 0.65(A:^//o)^^^cm^s) dimensions 
[T7] . Dimensionless time, to = lot^ and space, xq = x(/o/^)^^^, variables were introduced. For diffusion in water 
(/c = 10~^m^s~^) and nominal values ks = 10~^^cm^s~^ and Q = 10^^ cm~^s~^, we obtain Iq = ^/2 x 10~^s~^, 
simulating a 2d {3d) shallow (deep) pool with linear dimensions of / ~ 200 (50) cm. For the purpose of illustration, 



explicit results quoted below in Section III were computed using these values. 
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FIG. 2: Evolution of 2d chiral domains. Red (+1 on the color bar) corresponds to the L-phase and blue (-1 on the color bar) 
corresponds to the D-phase. Time runs from left to right and top to bottom. Top left, the near-racemic initial conditions. Top 
mid and top right, evolution of the two percolating chiral domains separated by a thin domain wall. Bottom left, environmental 
effects break the stability of the domain wall network. Bottom right, subsequent surface-tension driven evolution leads to a 
enantiomerically-pure world. 




FIG. 3: Time evolution of the spatially- averaged net chirality corresponding to the snapshots shown in figure [2] Stars denote 
times for snapshots and vertical lines mark the beginning and end of stochastic environmental influence. 
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FIG. 4: Punctuated Chirality. Impact of environmental effects of varying duration and fixed magnitude {a^ /al = 0.96) on the 
evolution of prebiotic chirality in 2d. Short events (last from left), which have little to no effect, should be contrasted with 
longer ones, which can drive the chirality towards purity and/or reverse its trend. (See, e.g. the green line.) 



1. Numerical Implementation 



The numerical implementation used a finite-difference leapfrog method with temporal step 6t = 0.005 and spatial 
step Sx = 0.2. Periodic boundary conditions were adopted on a 1024^ (256^) squared (cubic) lattice in 2d {3d). 
Near-racemic initial conditions were prepared in two steps. First, with S =1, the equation of motion for the net 
chirality A{t,x,y, z), eq. [2j was solved with a simple quadratic potential, V{A) = (1/2)^^ and low noise (a^ = 0.1) 
until the volumed- averaged chirality {A{t) = (l/V) J A{t^ x^y, z)dV) achieved a steady-state Gaussian distribution 
centered in the racemic state. By steady-state we mean that d{A{t))/dt ^ 0. Then, the noise {w{t^x^y^ z)) was 
turned off and the potential was switched to that of eq. [2l A{t^ x, y^ z) relaxed deterministically toward its dynamical 
equilibrium state, forming domains as illustrated in Fig.^ Of all such configurations, only those allowing for domain 
coexistence with near-racemic values were selected, about 38% of our total sample. That is, we discarded all cases 
where A{t) ^ ±1 quickly, as these would not be sensitive to low noise environmental disturbances. The final 100 
ensemble elements behaved typically as illustrated in Fig. [4] a coexisting domain structure was formed, achieving a 
near-static configuration for net chiralities not too far from the racemic state. At this point, the noise was turned on 
(at t > 600 in Fig. |4| in the equation for A{t^ x, ?/), simulating the stochastic impact of environmental influences with 
different intensity and duration. 



III. RESULTS 



In Fig. [4| we show several 2d runs where the environmental effects vary in duration, while their magnitude was 
set at a^/o^ = 0.96. Each colored line represents a prebiotic scenario, with sequentially occurring environmental 
disturbances of different duration. In order to investigate the impact of environmental effects on chiral selectivity, the 
scenarios reflect situations where there is no chiral selection, that is, where the two phases may coexist in dynamical 
equilibrium (mathematically, when d{A{t))/dt for A{t) ^ ±1; chemically, in a steady state). We observe that long 
disturbances can drive the net chirality towards purity ((^(t)) ^ ±1 for large t). Furthermore, note that subsequent 
events may erase any previous chiral bias, favoring the opposite handedness. In other words, environmental effects of 
sufficient intensity and duration can reset the chiral bias. This is true even if the system evolved toward homochirality 
prior to any environmental event. 

In Fig. |5j we summarize the results of a detailed statistical analysis of 100 2d runs that led to initial domain 
coexistence, that is, d{A)/dt ^ (see Fig. |4] for t < 600). The horizontal axis displays the magnitude of the 
disturbance in units of a^. The vertical axis gives the fraction of homochiral worlds, that is, those that after the 
disturbance obtain chiral purity. The colors represent the duration of the event. For > 0.96a^, that is, near the 
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FIG. 5: Fraction of 2d homochiral worlds after a single environmental event. Colors indicate duration of the event, while the 
horizontal axis labels its magnitude in units of a^, the critical value for chiral symmetry restoration. The statistical sample 
included 100 runs and the time scales of the events are in units t(/^^). 



critical region, all but the shortest events {t < 50/^^ ~ 1.5 months, for the nominal value of Iq = a/2 x 10~^s~^ 
mentioned previously) lead to statistically significant chiral biasing. Results in 3d are qualitatively very similar, 
although due to heavy CPU demand we limited our analysis to 50 short runs. 

i.e. those that 



It is interesting to note that the cases discarded from our statistical studies (see section II C 1 



spontaneously evolve toward homochirality from near-racemic initial conditions {A{t) ^ ±1 quickly) represent 62% 
of the total data set. Comparing to the fraction of systems evolving toward A{t) ^ ±1 after an environmental 
disturbance (see Fig. [5|, we see that long and high-intensity events are more efficient at evolving the system toward 
homochirality than systems starting from an initially random near-racemic state. In other words, large disturbances 
can not only drive a near-racemic system toward homochirality but also can reverse the chirality of homochiral 
solutions, effectively erasing any previous chiral signature (see, e.g., green line in Fig. [4|. 

These results suggest, in the one extreme, that the early Earth may have played host to numerous abiogenetic 
events, only one of which ultimately led to the Last Universal Common Ancestor through the usual processes of 
Darwinian evolution. This is consistent with work indicating that life may have become globally extinct more than 
once [? ]. In the other extreme, one may consider, at the very least, that biological precursors certainly interacted 
with the primordial environment and may have had their chirality reset multiple times before homochiral life first 
evolved. In this case, our results show that separate domains of molecular assemblies with randomly set chirality may 
have reacted in different ways to environmental disturbances. A final. Earth-wide homochiral prebiotic chemistry 
would have been the result of multiple interactions between neighboring chiral domains under mechanisms already 
described elsewhere |5l [TB] [18] . 

It is interesting to compare this scenario with the suggestion that life originated in hydrothermal vents [? ]. One 
could argue that spatially separate vents consist of independent ecosystems, each thus capable of evolving life with 
its own chirality. And yet, what is seen today is a dominance of the usual situation with L- amino acids and D- 
sugars. According to our results, if life originated in distinct hydrothermal vents, it would indeed have been possible 
for its precursors to have evolved into distinct chirality. One would then have life of different chirality in different 
hydrothermal vents. Does this offer a counterexample to our scenario? Not at all. Firstly, there are several obstacles 
for the prebiotic synthesis of organic compounds or polymers in hydrothermal vents [? ]. Second, even granting the 
unlikely possibility that life could have originated several times over in distinct hydrothermal vents where temperatures 
reach 350^ C, one must keep in mind that the whole ocean passes through them every ten million years. What happened 
in the distant past does not stay frozen in time. There is extensive mixing and recombination, leading to the exposure 
of separate chiral domains to each other and thus to competition and final preponderance of one handedness through 
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a prebiotic version of Darwinian evolution through natural selection. 



IV. CONCLUSIONS 

The possibility that environmental effects determined the homochirality of terrestrial life leads to several significant 
consequences. Any initial bias from extraterrestrial CPL, even if effective in the short term, would have been erased 
by subsequent environmental disturbances. Our results suggest that the detailed evolution of prebiotic chiral reaction 
networks is enmeshed with Earth's environmental history. Prebiotic Earth was a highly unstable medium, subject 
to active volcanism and meteoritic bombardment. Although the present analysis must be extended to describe such 
events in more quantitative detail, it is clear that prebiotic chirality, just as with the evolutionary history of living 
organisms, responded to environmental changes, alternating periods of stasis and violent activity. Our analysis predicts 
that other planetary platforms in this solar system and elsewhere could have developed an opposite chiral bias. As a 
consequence, a statistically large sampling of extraterrestrial stereochemistry would be necessarily racemic on average. 
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